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Topological materials and their unusual transport properties are now at the 
focus of modern experimental and theoretical research. Their topological properties 
arise from the bandstructure determined by the atomic composition of a material 
and as such are difficult to tune and naturally restricted to < 3 dimensions. 

Here we demonstrate that n-terminal Josephson junctions with conventional 
superconductors may provide a straightforward realization of tunable topological 
materials in n — 1 dimensions. For n > 4, the Andreev subgap spectrum of the 
junction can accommodate Weyl singularities in the space of the n — 1 independent 
superconducting phases, which play the role of bandstructure quasimomenta. 

The presence of these Weyl singularities enables topological transitions that are 
manifested experimentally as changes of the quantized transconductance between 
two voltage-biased leads, the quantization unit being 2e^/(7r/i). 

*To whom correspondence should be addressed; E-mail: roman-pascal.riwar@cea.fr 

Josephson junctions created by coupling two superconductors through a weak link have 
been studied extensively for many years^^-. The current across a Josephson junction yields 
information about the Andreev bound states (ABS) forming at the junction^"-. In turn, the 
ABS spectrum is determined by the properties of the junction and the superconducting leads. 
For instance, if the leads are topologically non-trivial, the Josephson effect may be used to 
probe these topological properties. In particular, the 47r-periodicity of the supercurrent 
indicates the presence of topologically protected zero-energy Majorana states^i^— , which 
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may arise in ID spinless p-wave superconductors, semiconductor nanowires with proximity- 
induced superconductivity, or at the surface of bulk materials. 

In this paper, we show that mr</t«-terminal Josephson junctions may be topologically 
nontrivial even if the superconducting leads are topologically trivial and no exotic materials 
are used to make the junction. Thus, the junction itself may be regarded as an artihcial 
topological material, which displays Weyl singularities. 

The rapidly growing field of 3D Weyl semimetals deals with a bandstructure that exhibits 
conical energy gap closings: Weyl points^^— . Unlike the Dirac point in graphene^^ that 
may be gapped out through an appropriate coupling, isolated Weyl points are topologically 
protected. They can be regarded as monopoles with ± charge. A topological invariant - 
the Chern number - defined on a surface in momentum space characterizes the total charge 
of the monopoles it encloses. These monopoles give rise to many unusual features, such as 
chiral edge states and associated surface Fermi arcs^^. Topological protection guarantees 
that the only way to induce a gap is either to annihilate two Weyl points of opposite charge 
by bringing them together, or to couple two cones at a finite distance in momentum space 
through breaking of momentum conservatioi>i^. 

Here, we demonstrate Weyl singularities in multi-terminal junctions with conventional 
s-wave superconducting leads, namely when the energy of the lowest ABS goes to zero at 
certain values of the superconducting phases such that the gap in the spectrum closes. Below, 
we also show that their topological property can be easily probed by the transconductance 
between two voltage-biased leads, which is proportional to the Chern number. 

We consider a junction with n superconducting leads connected through a scattering 
region (see Fig. [IK). The leads a = (0,1,..., n — 1} have the same gap A, though they may 
differ in the phase of the superconducting order parameter, (/>„. Due to gauge invariance, 
only n — 1 phases are independent, hence we may set 0o = 0. Likewise we choose to focus 
on a short scattering region, characterized by an energy- (and spin-) independent scattering 
matrix S in the basis of A = transport channels, Na being the number of channels 

in contact a. We also assume time-reversal symmetry, such that = S. 

The junction hosts a set of spin-degenerate Andreev bound states, indexed by k, whose 
energies > 0 are determined from the equation^ 


det 


1 - 


= 0 , 


( 1 ) 
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where x = arccos(£^/A), and is a diagonal matrix that assigns to each channel the phase 
factor of the corresponding terminal. The spin-degenerate ABS energies £’^(0) (see Fig. Hb) 
are periodic in all phases with a period 27r. The total phase-dependent energy of the jnnction 
reads E = — l/2)£fc, Uka = 0,1 being the occupation of the state k with spin a. 

If the gap closes at certain values of the phases, the spectrum is doubly degenerate 
at £ = 0. We can then describe the lowest energy band in the vicinity of the gap closing by 
the two-by-two Weyl Hamiltonian (Supplementary Discnssion), 


£w ^ ^ hjTj , hi ^ ^ j (2) 

i=x,y,z a 

with the Pauli matrices in the basis of the two degenerate states corresponding to 

eigenvalne £ = 0. The fields hi depend linearly on (50 = 0 — through the real matrix 

M. 


The form of the Weyl Hamiltonian in Eq. (EJ indicates that we need at least three 
parameters to tune the system to the degeneracy point, hj = 0. Thus, for 4 terminals 
with 3 independent phases, the Weyl singnlarities appear as points in the 3D phase space. 
For 5 terminals, the Weyl singnlarities occur in general as ID cnrves in the 4-dimensional 
space of phases. This opens np the possibility to realize a topological material in arbitrary 
dimensions. Note that such multi-terminal Josephson junctions cannot be characterized by 
means of the standard “periodic table” of topological semimetalsi^'i^, due to the distinct 
behaviour of the quasimomenta 0 under particle-hole symmetry. A proper classification 
could be envisioned along the lines of Ref.—. 

The topology of the jnnction can be characterized by a set of Chern numbers. A Chern 
number may be defined in the 2D snbspace of two phases (f)a and 0/3, through the local Berry 
curvature of the ABS. The Berry curvature for the bound state k with spin a is related to 
its wave function \(pka) through 


Bf = -2Im 


d^k<T d(pka \ 
90a 90/3 / ■ 


( 3 ) 


Note that does not depend on spin (Supplementary Discussion). The total Berry 
cnrvature of the many-body superconducting state then reads £“^ = ’Yhkai'^ka — V2)£fc^; 
in analogy to the expression for the energy. For fixed occupations = 0,1, the integral 
of the Berry curvature over the elementary cell yields an integer, the Chern number = 
£“^/(27r). 
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FIG. 1. General setup of the multi-terminal junction, and examples of typical ABS 
spectra, a) The superconducting leads with phases a = 0,..., n — 1, are connected through 
a scattering region described by the scattering matrix S. b) Generic ABS energy spectrum versus 
(pi, away from a Weyl singularity, c) ABS energy spectrum versus cpi , where the other phases are 
tuned to a Weyl singularity. Note the gap closing. 

Since the Weyl singularities appear as points in the 3D phase space, a third phase (p^ 
may be used to tune the system through Weyl points, thus changing the Chern number. We 
see that a given state k only contributes to when it has an even occupation number, 
The topological signature is therefore robust provided the fermion parity is 
preserved, a fact which has also been pointed out for other topological systems realized in 
superconductors^. Quasiparticle poisoning, by which a single quasiparticle enters a bound 
state, enables transitions between even and odd ABS occupatiorr^i, and can therefore affect 
the signal. In the following, we focus on the ground state, rika = 0 for all k and a, where 
the Chern number may be nontrivial. 

Importantly, the current response of the junction with slowly varying phases reveals the 
Chern number. Biasing lead (3 with a voltage eV^ <C A gives rise to the instant current to 
contact a (Supplementary Discussion) 

where = 2eVj3/h. The first term corresponds to the adiabatic current and the second 
term is the first order correction in the phase velocity. Let us now apply constant voltages 
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to two leads. For incommensurate voltages, the two phases uniformly sweep the elementary 
cell. In the dc limit, the adiabatic current contribution then averages out, and the Berry 
curvature is replaced by its average value. Thus, we find that the dc current is linear in the 
voltages, and the transconductance is defined by the Chern number 

2p2 

with =-(5) 

irh 

Equation ([SD shows that multi-terminal junctions exhibit a dc current response typical for the 
quantum Hall effect, although based on different physics. The transconductance quantum is 
four times bigger than in the quantum Hall effect, which can be traced to the 2e charge of the 
superconducting Cooper pairs. To extract the small dc signal, the averaging time needs to 
be sufficiently long. The relevant time scale is determined by the low-frequency current noise 
(Supplementary Discussion). Note that quantum Hall-like conductance quantization has also 
been proposed in superconducting devices with finite charging energy and hosting quantum 
phase slips^^. Furthermore, superconducting junctions with a gate-tunable charging energy 
may realize topologically-protected discrete charge pumping^S. 

We now focus on a 4-terminal junction and investigate the energy spectrum as a function 
of the three independent phases 0 p 2 , 3 - As mentioned above, such a 3 D bandstructure may 
host Weyl points with positive or negative topological charge. The Nielsen-Ninomiya theo- 
rem 2 i implies that the total topological charge of the system is zero, such that the number 
of Weyl points is always even. Furthermore, time-reversal invariance corresponds to a map¬ 
ping from (j) to —0, hence a Weyl point at has a counterpart at —with the same 
topological charge. Thus, Weyl points come in groups of 4. 

In the simplest case, where each contact contains only one channel, the system may 
realize 0 or 4 Weyl points, corresponding to a topologically trivial or nontrivial 3D material, 
respectively. If a scattering matrix yielding Weyl points is found, small changes in the 
scattering matrix only modify their position, but cannot gap them. Namely, as the Weyl 
points carry a topological charge, individual Weyl points are stable and annihilation is 
possible only when two Weyl points with opposite charges coincide. 

A specific example is shown in Fig. [21 The position and charge of the 4 Weyl points 
is shown in Fig. [ 2 ^ 1 . Without loss of generality, we fix the phase 03 and compute the 
transconductance between voltage-biased contacts 1 and 2. In Fig. [2J3, one can clearly 
see that the transconductance increases (decreases) by 2 e^/ (nh) when 03 passes through a 
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FIG. 2 . Topological characterization of the 4-terminal junction for the single-channel 
case, a) Position of the four Weyl points in the space of 4>i^2,3 of the single-channel 4-terminal 
junction, the colour code indicating the respective charge, b) The resulting transconductance 
indicating the Chern number, as a function of ^3 for the same single channel junction as in panel 
a). 

Weyl point with positive (negative) charge. Interpreting ^3 as a control parameter rather 
than a quasimomentum, we thus see that the 2D bandstructure of the system as a function 
of (pi and (p 2 undergoes a topological transition when (p^ passes through a Weyl point. 
The transconductance directly measures the Chern number characterizing the corresponding 
2 D topological phase. Note that the transconductance satisfies the relation ^ 3 ) = 

—G^^{(p 3 ) due to time-reversal symmetry. 

By randomly generating scattering matrices from the circular orthogonal ensemble^^, we 
find that about 5% of scattering matrices give rise to 4 Weyl points (Supplementary Dis¬ 
cussion). More Weyl points can be obtained in multi-channel junctions where the maximal 
number of Weyl points is roughly proportional to the number of channels, and the proba¬ 
bility to have no Weyl points is small. As a consequence, a greater variety of 2 D topological 
phases with higher Chern numbers can be realized in that case. This is shown in Fig. E] for a 
multi-channel junction hosting 36 Weyl points, where the maximal Chern number is 3. Re¬ 
cently realized few-channel cross junctions^^ are promising to observe the transconductance 
in a 4-terminal junction. 


(p2 
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FIG. 3. Topological characterization of the 4-terminal junction for the multi-channel 
case. Chern number as a function of ^3 for a multi-channel 4-terminal junction, where the contacts 
1,2,3, and 0 contain 12,11,10, and 9 channels, respectively. In this particular example, the junction 
hosts 36 Weyl points. 

We now turn to 5-terminal junctions. In that case, the Weyl singularities appear as 
closed loops in the 4D space of phases. The simplest way to visualize them is to consider 
the additional phase ^4 as a tuning parameter of the 3D system described by the phases 
( 01 ; 02 , 03)- Tuning 04 the Weyl points move, but remain at zero energy. Note that in the 
3D subspace, time-reversal symmetry is effectively broken, as for a fixed nonzero 04 , a Weyl 
point at (0i°\ 02 °\ 03°^) does not have a counterpart at — (0i°\ 02 °\ 03 °^) anymore. The only 
constraint at a finite 04 is that the number of Weyl points is even. Once two Weyl points 
with opposite charge meet, they annihilate. Thus their trajectories describe closed loops in 
the 4D space of all phases. 

As before, the presence of Weyl singularities may be probed by the transconductance 
between two voltage-biased terminals, say terminal 1 and 2 . As a function of the two other 
phases 03 and 04 , we now obtain areas with different quantized values of the transcon¬ 
ductance, corresponding to different Chern numbers. The boundaries of these areas are 
given by the projections of the Weyl loops on the (03, 04)-plane. An example is shown in 
Fig. m Here time-reversal symmetry of the scattering matrix manifests itself in the relation 
G^^( — 03 ,—04) = —G*^^(03, 04)- 

The above considerations are straightforwardly extended to a larger number of terminals. 
Furthermore, as we saw for the 5-terminal junction, breaking time-reversal invariance does 
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FIG. 4. Topological characterization of a 5-terminal junction, each contact having a 
single channel. The colored areas display a nonzero Chern nnmber The bonndaries of these 
areas correspond to the projection of the Weyl singularity lines to the (<^ 3 , <^ 4 )-plane. In a 3D 
subspace the curves can be assigned a charge. 


not lift the topological protection of the Weyl points at zero energy, as long as spin degener¬ 
acy is preserved. When spin degeneracy is lifted due to a Zeeman field or due to spin-orbit 
interactions, we expect the Weyl points to shift away from zero energy while remaining 
stable. The possibility of realizing zero-energy states in 3-terminal junctions with strong 
spin-orbit interaction was studied in Ref.—. 

To summarize, we predict the existence of topological Weyl singularities in the ABS 
spectrum of n-terminal superconducting junctions with n > 4. These Weyl singularities 
manifest themselves in a quantized transconductance between two voltage-biased contacts, 
when the remaining phases are tuned away from the singularities. Tuning the system through 
a Weyl singularity, the conductance displays a step, signaling a topological transition of 
the 2 D system described by the two phases of the voltage-biased contacts. As with all 
topological properties, this signature is robust provided the system remains in its ground 
state. Thus, the quantized transconductance should be easily accessible experimentally at 
low temperatures and voltages. 
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I. DERIVATION OF THE WEYL HAMILTONIAN 


In the following, we provide the derivation of the Weyl Hamiltonian, Eq. (2) of the main text, 
and associated properties of the Andreev bound state (ABS) spectrum. 

The topological properties of multi-terminal junctions result from the appearance of Weyl sin¬ 
gularities in the ABS spectrum. The ABS for an n-terminal superconducting junction are deter¬ 
mined through^ 



( 

S 

0 


( •m\ 

0 e'^ 

, with Sn = 

O 

5% 

and SA = 

O 

T 

_^ 




(SI) 


Here, the two-by-two structure indicates the Nambu space of electrons and holes. The scattering 
matrix S ^ describes the normal metal scattering region, where S and S * provide the scattering 
amplitudes for electrons and holes, respectively. The matrix S a accounts for the phase acquired in 
the Andreev reflection processes, where = arccos(£'/A) and is a diagonal matrix that assigns 
to each channel the phase of the corresponding terminal. Note that the wave function \i}/) is written 
in terms of the outgoing states in a given spin sector cr =t, f. 

To determine the ABS energy bands, it is sufficient to reduce the above equation to the electron 
subspace, where we find the determinant condition for the bound state energies E as 


det 1^1 - e = 0 , 


(S2) 
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with the unitary matrix A(0) = Se“^S*e~“^. The matrix A possesses the particle-hole symmetry 
A* = U^AU, where U(^) = . This implies that the eigenvalues of A come in pairs^, 

corresponding to energies +Ek, where = Acos(aicl2) with 0 < < n. We assign eigenvectors 

I'Fp and |'P^) = C/|'I'P* to the pair of eigenvalues respectively. 

If there is a zero-energy solution at the spectrum of the unitary matrix has a dou¬ 

bly degenerate eigenvalue -1. The corresponding orthogonal eigenvectors are given as \a^) = 
and \a~) = In the vicinity of the singularity at we expand the determinant 

equation (IS2h for small 50 = 0 - «: 1 and E <s: A. Up to first order, we find « -1 + 2iEI A 

and 

A(0) « A(0^°^) + iS 6^S . (S3) 

Projecting Eq. (IS2h onto the subspace defined by \a^) and \a~), and keeping only the lowest-order 
terms, we find the determinant equation 


det [E - //w] = 0 . 


(S4) 


It defines the eigenvalue problem for the lowest band, described by the two-by-two Weyl Hamil¬ 
tonian, 


(a+l 50 - 5505+ |a+) (a+| 50 - 5505+ la") 

(a-|50 - 5505+|a+) (a'l 50 - 5505+la") ^ ' 

Using the particle-hole symmetry \a~) = Eq. (IS5]) may be simplified to 



(S5) 


^ , (S6) 

i=x,y,z 

where Tx,y,z are Pauli matrices in the basis of da"^), |a“)}, and the real prefactors take the form 
hx + ihy = A<50)_+ and = f ((50)++ - <50)__). As a consequence, the energy spectrum is linear 
in 50 and reads E = ± + hy + h\. We see that 3 independent parameters are necessary to 

tune the energy to zero. Thus, the Weyl singularities in general appear as points in the space of 
three independent phases of the 4-terminal junction, and as curves in the space of four independent 
phases of the 5-terminal junction. 

In a given 3D subspace, one may assign a charge to each Weyl point. The sign of a Weyl 
charge is most conveniently defined by rewriting the Weyl Hamiltonian in matrix form, /fw = 
Yjai S(f>aMaiTi with Mq,, = ds,p„hi. Then, in a subspace of three independent phases, M is a square 
matrix, and the sign of the charge is given by the sign of the determinant, det[MQ,,]. 
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While particle-hole symmetry is local in time-reversal symmetry links ^ and As a 
consequence, a junction with time-reversal symmetry gives rise to pairs of Weyl singularities at 
and Indeed, using = 5, we find the relation A(-0) = e~‘^A\^)e‘K Using this 

relation, we find that the Weyl Hamiltonians near and take the same form and hence, in 
a given 3D subspace, the corresponding Weyl points have the same charge. 

We note that there is the possibility of Weyl singularities at finite energy. Due to particle-hole 
symmetry, they appear in pairs at energies ±E with the same charge. Because of particle-hole 
symmetry and spin degeneracy, solutions at energy -E can be ascribed to solutions at energy E 
in the opposite spin sector—. Thus, zero-energy Weyl points are singly degenerate, while finite 
energy Weyl points are doubly degenerate. As the finite energy Weyl singularities do not affect 
the ground state Chern numbers used to characterize the system, we do not discuss them any 
further. Note that the outgoing wavefunctions in particle-hole space used to compute the Chem 
numbers are given as Furthermore, \(//ko-) is identified with in the 

corresponding spin sector cr. 


II. DERIVATION OF THE CURRENT 


In the following, we establish the connection between the current and the Berry curvature for 
phases that change slowly in time. 

The current operator through lead a is defined as 

dH 

la = 2e— , (S7) 

0(pa 

where H is the Bogoliubov-de Gennes Hamiltonian. In order to calculate its expectation value 
for time-dependent phases 0(0, we introduce the basis of instantaneous wave functions of the 
time-dependent Bogoliubov-de Gennes Hamiltonian H(t), such that Ek(t)\(pkcr(0) = H(t)\(pka-(0)- 
Solving the time-dependent equation ifi\(p) = H(t)\(p) in that basis, up to first order in phase velocity 
0, we obtain the current contribution from state k with spin cr as 


lakit) « 2e 


1 dEkit) 
% d(pa 


,d {^^^(^(Ol 

50a 


|0/:(t(O) "I" ^X0/:(t(OI 


5 |^A;(r(0) 
d(l>a 


(S8) 


Note that lak does not depend on spin. The first term corresponds to the adiabatic supercurrent. In¬ 
troducing the Berry curvature = -2Im ^d^^{cpko-\d^^\(pkcr)\, the second term reads -2e 2/j 0/3^*^- 
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The many-body expectation value of the current is computed as 


kcr 



(S9) 


where Uko- = 0,1 is the occupation of the ABS k with spin cr. 

Equations (IS8l) and (IS9l) can be used to compute the adiabatic correction to the current, if 
the energy spectrum is discrete and the occupations of each state do not change with time. For 
a few-channel junction, the level spacing between the ABS scales as ~ A, thus providing the 
adiabaticity condition A. On the other hand, defining the current contribution of the states 

in the continuum above the gap, as well as of ABS that reach the gap edge at some values of 
the phases, is problematic. In particular, adiabaticity would be broken at any phase velocity for 
the states in the continuum, while ABS that reach the gap edge may lead to stricter adiabaticity 
conditions (see next section). In the following, we argue that those states do not contribute to the 
topological properties considered in our work. 

Indeed, multi-terminal junctions are known to be topologically trivial in certain cases. For 
instance, let us consider the case where terminals a and jS are voltage biased. When all the other 
phases (py (y 4^ a,/5) are set to zero, the junction is effectively a 3-terminal junction which is 
topologically trivial. In that case, the associated Chern number vanishes. A change of C“^, 
corresponding to a topological transition, may occur when a state crosses the Fermi level as the 
phases py are varied. As the energy spectrum is gapped, this can only occur for the lowest energy 
states. Thus, only those states contribute to the topological properties, and in the following we 
only retain their contribution in Eqs. (IS^ and (IS9l) . This establishes the regime of validity for 
the relation between the transconductance and Chern number given in the main text, see Eq. (5). 
To compute the Chern numbers, we use the outgoing wavefunctions lipko-) defined in the previous 
sectional. 


III. COMMENT ON ANDREEV BOUND STATES REACHING THE GAP EDGE 

Fet us return to the possibility of the ABS with the highest energy reaching the gap edge. 
Indeed we find that, for a system with an even number of total channels, this occurs along (n- 1)- 
dimensional subspaces in the space of n independent phases^. 

An argument for this fact can be given in analogy with the derivation of the Weyl Hamiltonian. 
Namely, an ABS reaching the gap edge at some (^g requires that the unitary matrix A(^g) has a 
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FIG. SI. Zoom on the highest ABS of the spectra shown in Figs, lb and Ic in the main text. The gap 

edge touchings are clearly visible. 


doubly degenerate eigenvalue +1. We denote the corresponding orthogonal eigenvectors \b'^) and 
\b~) = U(^g)\b'^y. Expanding Eq. (IS2h for small 6^ = ^ - 0g, we can reduce the determinant 
equation to 


det 


signCE) 


8(A-|E|) 


+ 


(<d^)++ - (d0)—) 


- T, 


= 0 


(SIO) 


where is a Pauli matrix in the basis {|^^), \b )}. 
Thus, the solution for the ABS energies reads 


E 


+A 




(Sll) 


Furthermore, we see that ABS wave function has a discontinuity at 0g. Namely, the wave function 
of the state at positive energy is given by \b'^) for {6^)++ - (6^)— < 0 and by \b~) for (5^)++ - 
{(50)__ > 0, and vice versa for the states at negative energies. 

In contrast to the discussion for the zero-energy Weyl points, here a single parameter is sulR- 
cient to tune the ABS energy to the gap edge. Therefore, the highest ABS reaches the gap edge 
along (n - l)-dimensional subspaces in the space of n independent phases. This is visible in 
Figs, lb and Ic in the main text. A zoom is shown in Fig.|SU 

When sweeping two phases, the curves where the highest ABS reaches the continuum can¬ 
not be avoided. If this state was occupied, this would break adiabaticity even for arbitrary slow 
driving. (Introducing a finite energy-dependence of the scattering matrix or spin-orbit interactions 
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may create a finite level repulsion between the highest ABS and the continuum and, thus, restore 
adiabaticity, possibly under stricter conditions. This would need further investigation.) However, 
as far as the topological characterization of the junction is concerned, we came to the conclusion 
that this peculiarity should not change our results. 


IV. ESTIMATION OF THE CURRENT NOISE 


The topological signal of interest is a small dc contribution on top of a large ac adiabatic current. 
This makes it important to provide an estimate for the measurement time Tq required to sufficiently 
average the current signal. In particular, the lower bound for Tq is provided by S //|7p, where 7 
denotes the dc current signal and S/ is the low-frequency current noise. The dominant contribution 
to 5/ stems from fluctuations of the phases due to voltage noise and will be estimated in the 
following. 

We define the current noise as 


5/ = 


f 


dr + r)) , 


(S12) 


where the brackets (...) denote the expectation value with respect to the statistical ensemble. 
As the adiabatic supercurrent, i.e., the first term in Eq. (IS^ . is much larger than the correction 
proportional to the Berry phase, it dominates the noise. We exploit the fact that this current is 
periodic in both phases 4>a and (pp, such that it may be expanded in Fourier harmonics. 


I(<Pa,<Pls)= 

na,np=-oo 

Substituting Eq. (IS13h into Eq. (IS12h . we obtain 


I dT E I/.. |2^^ma[0a(r)-0a(r+T)] ^ 

na,n/3 


(S13) 


(S14) 


as the two phases fluctuate independently and the ac contribution averages to zero. 

In order to evaluate the averages, we assume white noise in the voltage sources, described by 




(S15) 


with the rate Eq, associated to the voltage noise in lead a. Using (pa{t + t) - (paiO = + 

(lefK) ds6Va(s), where oja = leValK we find 


{^gina[<Pa{t)-(pa{t+Ty\\^ _ ^ 


-inaa>ar-Tanl\T\ 


(S16) 
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and subsequently 


SI — 2 


"a^a + T; 




2\2 




(UaiOa + npiOpY + (r„n2 + Yi^np 


(S17) 


We see that S / depends on frequencies co ~ coaj} in a quite complex fashion. In the limit of 
low voltage noise, » F with F ~ ^a,p, we have to distinguish between commensurate and 
incommensurate voltages Vaj}. In the commensurate case, the first term in the denominator of 
Eq. (IS17h vanishes for a pair {Ua, Uj^), yielding 5/ oc 1/F. By contrast, in the incommensurate case, 
the first term in the denominator of Eq. (IS 171) dominates, yielding a much weaker noise, S/ oc F jaP-. 
Thus, the current noise strongly varies with co. Averaging over a window of width ~ oj yields the 
average noise Sjocl jco. On the other hand, in the noisy regime cu ^ F, we find that 5/ oc 1/r, 
irrespective of the voltages being commensurate or not. Thus, the current noise is actually weaker 
than the average current noise in the opposite regime. Furthermore, incommensurability is no 
longer required for averaging, as the strong phase fluctuations take care of covering the entire unit 
cell. 

This brings us to the conclusion that the strong voltage noise regime ru <?c F may be favourable 
for averaging - provided, of course, that the noise is still in the limit %Y A. To estimate the 
measurement time, we use ~ and I ~ e^V/fi with V ~ Vajs, together with the estimate 
for SI found above, to obtain 


To> 


A^ 1 
(eVp f ■ 


(S18) 


V. OCCURENCE OE WEYL POINTS EOR 4-TERMINAL JUNCTIONS 

Here, we estimate how often random scattering matrices give rise to Weyl points at zero energy. 

As discussed in the first section of the supplementary, such Weyl singularities emerge when 
the matrix A has the eigenvalue -1. Whether this occurs indeed depends on the properties of the 
junction, encoded in the scattering matrix S. We performed a numerical analysis for 4-terminal 
junctions, both in the single-channel and multi-channel cases. For this purpose, we generated 
symmetric scattering matrices according to the circular orthogonal ensemble using the following 
recipe^: We randomly generated hermitian matrices TF from a Gaussian ensemble, and numer¬ 
ically diagonalized them, 7F = where FF is a unitary matrix and D is a real diagonal 

matrix. A random symmetric scattering matrix was then generated as 5 = For each S , we 

numerically checked the existence of zeros of the function |det[l + A(0)]| in the space of phases. 
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FIG. S2. Histogram displaying the occurrence of random scattering matrices yielding N Weyl points 
for the four-terminal multi-channel junction. 

In the single-channel case, we ran the check for 965 randomly generated matrices S, out of 
which 46 gave rise to zero-energy solutions. Thus, we found that a total of roughly 5% of all 
scattering matrices yield Weyl points, while the remaining 95% provide a trivial junction. 

When increasing the number of channels in each terminal, we found that the maximal number 
of Weyl points scales with the number of channels, and that the probability of a junction without 
Weyl points decreases significantly. A total of 324 random scattering matrices were generated for 
a junction with four terminals, where the terminals have 12,11,10, and 9 channels, respectively. 
In Fig. 1^ we show the histogram displaying the occurrence of randomly generated scattering 
matrices that provide N zero-energy Weyl points. Only 4 scattering matrices gave rise to a junction 
without zero-energy Weyl points. Note that our algorithm has a small, but finite probability to miss 
some zeros. As a consequence, while the number of Weyl points has to be a multiple of 4 because 
of time-reversal symmetry (see main text), the algorithm misidentifies a few cases as having an 
odd multiple of 2 Weyl points. 


SI C. W. J. Beenakker, Phys. Rev. Lett. 67, 3836 (1991). 

For an odd number of channels, there is an additional eigenvalue 1 corresponding to a state at the gap 
edge. 

Note that in the main text we restrict ourselves to positive energies while keeping both spin sectors. 
Note that the Chern number may be computed using either the outgoing or the incoming wave functions 
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as they are related through Sn which does not depend on the superconducting phases. 

As pointed out previously, in the case of an odd number of total channels, there is a state at A for any 
value of the phases. 

M. L. Mehta, Random matrices, vol. 142 (Elsevier/Academic Press, Amsterdam, 2004). 


